A process-based, climate-sensitive model to derive methane emissions 
from natural wetlands: application to 5 wetland sites, 
sensitivity to model parameters and climate 


Bernadette P. Walter 1 * and Martin Heimann 2 


'Max-Planck-Institut fiir Meteorologie, Bundesstrasse 55, 20146 Hamburg, Germany (now at 
Columbia University/NASA Goddard Institute for Space Studies, 2880 Broadway, New York, NY 
10025, USA) 

2 Max-Planck-Institut fiir Biogeochemie, Postfach 100164, 07745 Jena, Germany 


submitted to: Global Biogeochemical Cycles, July 1999 


(*corresponding author: email: bwalter@giss.nasa.gov) 



Abstract 

Methane emissions from natural wetlands constitutes the largest methane source at present and 
depends highly on the climate. In order to investigate the response of methane emissions from 
natural wetlands to climate variations, a 1 -dimensional process-based climate-sensitive model 
to derive methane emissions from natural wetlands is developed. In the model the processes 
leading to methane emission are simulated within a 1 -dimensional soil column and the three 
different transport mechanisms diffusion, plant-mediated transport and ebullition are modeled 
explicitly. The model forcing consists of daily values of soil temperature, water table and Net 
Primary Productivity, and at permafrost sites the thaw depth is included. The methane model is 
tested using observational data obtained at 5 wetland sites located in North America, Europe 
and Central America, representing a large variety of environmental conditions. It can be shown 
that in most cases seasonal variations in methane emissions can be explained by the combined 
effect of changes in soil temperature and the position of the water table. Our results also show 
that a process-based approach is needed, because there is no simple relationship between these 
controlling factors and methane emissions that applies to a variety of wetland sites. The 
sensitivity of the model to the choice of key model parameters is tested and further sensitivity 
tests are performed to demonstrate how methane emissions from wetlands respond to climate 
variations. 

Introduction 

Methane is one of the important greenhouse gases, contributing about 22% to the greenhouse 
effect at present [Lelieveld et al., 1998]. Its atmospheric concentration has increased by a factor 
of 2.5 since the onset of industrialization and is now 1720 ppbv. Since methane has a large 
radiative effect - one unit mass of CH 4 has 2 1 times the radiative effect of one unit mass of 
C0 2 - changes in the atmospheric methane concentration affect the temperature on earth 
[IPCC, 1996]. Natural wetlands are the biggest methane source at present, contributing about 
40% to the total methane emissions and form the major non-anthropogenic methane source 
[Hein et al., 1997]. Furthermore, unlike most other methane sources, methane emissions from 
natural wetlands depend highly on the climate being influenced by temperature as well as by 
wetness. 

The emission of methane from natural wetlands is a result of biological and physical processes 
taking place in the soil: methane production by methanogenic bacteria under anaerobic 
conditions, methane oxidation by methanotrophic bacteria under aerobic conditions and 
transport of methane to the atmosphere [e.g. Conrad, 1989]. Numerous studies have been 
carried out to investigate which factors control these processes. It has been shown that the most 
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dominant determinants are: 


1. The position of the water table, because it determines the extent of the anoxic soil zone, 
where methane is produced, and the oxic soil zone, where methane is oxidized [e.g. Bartlett 
and Harriss, 1993; Bubier et al., 1995; Dise et al., 1993; Fowler et al., 1995; Funk et al., 
1994; Moore and ' Roulet, 1993; Morrissey and Livingston, 1992; Sebacher et al., 1986]. 

2. The soil temperature, since it influences the rates at which microbiological processes such 
as degradation of organic matter (i.e. the production of substrate for methanogenesis), meth- 
ane production and methane oxidation occur [e.g. Baker-Blocker et al., 1977; Bartlett and 
Harriss, 1993; Bubier et al., 1995; Christensen et al., 1995; Crill et al., 1988; Dise et al., 
1993; Fowler et al., 1995; Frolking and Crill, 1994; Kettunen and Kaitala, 1996; Morrissey 
and Livingston, 1992; Whalen and Reeburgh, 1992]. 

3. The availability and quality of suitable substrate for methane production. Several authors 
observed a correlation between factors indicating substrate availability, such as the Net Pri- 
mary Productivity (NPP), the Net Ecosystem Production (NEP) or the amount of biomass 
and methane production. Others observed enhanced methane emission after addition of sub- 
strate to the soil or from sites with a higher substrate quality [e.g. Bridgham and Richard- 
son, 1992; Klinger et al., 1994; Morrissey and Livingston, 1992; Valentine et al., 1994; 
Whiting et al., 1991; Whiting and Chanton, 1992; Whiting and Chanton, 1993], 

4. The pathways by which methane is transported to the atmosphere can be crucial for deter- 
mining the fraction of produced methane that is emitted into the atmosphere. Transport can 
proceed by molecular diffusion, ebullition and transport through the stems of vascular 
plants [Conrad, 1989]. Depending on the prevailing soil conditions the occurrence of one or 
more of those mechanisms can significantly alter the amount of emitted methane: 

• Ebullition can bring methane much faster to the water table than diffusion and thus increase 
methane emission [Bartlett et al., 1990; Boon and Sorrell, 1995; Devol et al., 1990; Wass- 
mann et al., 1992], 

• Plant-mediated transport can enhance methane emission through bypassing the often exist- 
ing oxic top soil layer and thereby avoiding methane oxidation there [Bartlett et al., 1992; 
Boon and Sorrell, 1995; Bubier et al., 1995; Chanton et al., 1992a; Chanton and Dacey, 
1991; Dacey and Klug, 1979; Happell et al., 1993; Holzapfel-Pschom et al., 1986; Morris- 
sey and Livingston, 1992; Schimel, 1995; Shannon et al., 1996; Shannon and White, 1994; 
Tom and Chapin, III, 1993; Waddington et al., 1996; Whiting and Chanton, 1992; Yavitt 
and Knapp, 1995]. On the other hand, atmospheric oxygen can be conducted by vascular 
plants down to the rooting zone. Therefore, in water saturated soils covered by vascular 
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plants a small oxic zone establishes around the root tips leading to methane consumption 
there [Gerard and Chanton, 1993; Holzapfel-Pschom et al., 1986; Schipper and Reddy, 
1996; Schiitz et al., 1989], 


Those controlling factors are not independent of each other, e.g. higher temperatures can lead 
to lower water table levels due to enhanced evapotranspiration. Thus, higher temperatures will 
lead only to higher methane emissions, if the soil stays water saturated. For example Roulet et 
al. [1992] found that temperature explained the variances in methane emissions as long as the 
soil was water saturated, but as soon as the water table decreased below the soil surface 
methane emissions were determined by the position of the water table. Valentine et al. [1994] 
observed that the temperature dependence of methane production increased with increasing 
amount and quality of substrate at sites that were substrate limited. Hence, since the 
environmental factors controlling the involved processes are not independent of each other and 
their relative importance in influencing methane emissions depends on the conditions present 
different factors are determining at different times and sites. Therefore, the results of these 
studies are site-specific and depend on the situations prevailing during the measurements. 

Some of the authors named above have also developed simple models to predict methane 
emissions from a wetland site. Those models are based on observed empirical relationships 
between the factors controlling methane emissions and methane fluxes and, as discussed 
above, are not capable of simulating methane emissions from other wetland sites and under 
conditions different from those prevailing when the data the models are based on were 
obtained. Besides the model described and used here there are two models in the literature 
using a more process-based approach, one by Cao et al. [1996] and one by Potter [1997], The 
model by Cao et al. [1996] simulates the carbon dynamics in the soil and derives the methane 
production rate as a function of the amount of decomposed organic carbon, the position of the 
water table and the temperature. Methane emission rates are calculated as the difference 
between methane production and oxidation rates. The model is applied on the global scale, but 
it has not been tested against data from any wetland site. The model by Potter [1997] is based 
on an ecosystem cycling model and simulates methane production rates from a microbial 
production ratio of C02-to-CH 4 that changes as a function of the water table depth. The three 
different transport mechanisms are modeled as simple functions of methane production, water 
table and ecosystem type. However, the model when tested against data from one field site 
does not reproduce observations very well. 

The objective of this work is to provide a model that can be applied to simulating methane 
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emissions from different wetland types in various regions as a function of the prevailing 
climate. The idea is to have a tool that can be used to study climate-dependent variations of 
methane emissions from natural wetlands on a larger (e.g. global scale) as it has already been 
done to investigate interannual variations of methane emissions from natural wetlands under 
present-day climatic conditions [Walter, 1998]. Therefore, the model used here is a more 
process-based, climate-sensitive model. In the remainder of this paper we give an extended 
description of a slightly modified version of the methane emission model briefly described in 
Walter et al. [1996] and show results of tests of the model against data from 5 different test 
sites located in North America, Europe and Central America, representing a large variety of 
environmental conditions. In addition, we present the results of sensitivity studies of the model 
and demonstrate how methane emissions from naturals wetlands might change under possibly 
changed climatic conditions in the future. 


Model Description 


Since the processes leading to methane emission from wetlands occur in the soil the following 
model structure is chosen (Fig. 1): a 1 -dimensional soil column is divided into 1 cm thick 
parallel layers. The position of the water table is assumed to be the boundary between the 
anoxic and the oxic soil zone. Methane is produced in the layers below the water table. The 
methane production rate is calculated as a function of the soil temperature and the Net Primary 
Productivity (NPP) which is assumed to be a measure for substrate availability. In the layers 
above the water table methane is consumed. The methane oxidation rate is calculated using the 
Michaelis-Menten equation. Besides, it is a function of the soil temperature. Three different 
transport mechanisms are modeled explicitly: (1) molecular diffusion through the water or air 
filled soil pores and the standing water, if the water table is above the soil surface, (2) transport 
by ebullition from the depth where bubbles are formed up to the water table and (3) transport 
through plants from soil layers above the rooting depth directly up to the atmosphere. The 
model forcing consists of daily values of the position of the water table, the soil temperature 
profile and the NPP. The model output are methane fluxes to the atmosphere and methane 
concentration profiles in the soil, both obtained on a daily basis by numerically solving the 1- 
dimensional continuity equation within the entire soil/water column: 


Figure 1 


| -CchS 1 > *) = + Qebull^’ *) + Qplanfr Z ) + R prod^ Z ) + R oxid^ Z ) 0 ) 

where the methane concentration at time t and depth z, F^t,z) the diffusive flux 
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of methane through the soil, Q e b u ii(t,z) and Q p i ant (t,z) represent sinks due to ebullition and 
plant-mediated transport, respectively. R pro( i(t,z) is the methane production rate, while 
Roxidft’Z) denotes the methane oxidation rate. In the following, each of the terms will be 
described in detail. 

Methane production rate R prod 

There are two factors controlling the rate of methanogenesis: (1) the availability and quality of 
suitable substrate which originates mainly from root exudates, dead fine roots and the input of 
plant litter and (2) temperature [e.g. Conrad, 1989; Heyer, 1990]. We consider only fresh 
organic matter produced by plants to be suitable substrate for methanogenesis [e.g. Chanton et 
al., 1995; Bridgham and Richardson, 1992] and assume the availability of substrate to be 
connected with Net Primary Productivity (NPP). Its variation with time f in (t) has been 
parameterized in the following way: it is calculated as a function of the variation of the NPP 
with tim e,f N pp(t), which is derived from the simulated NPP(t) of the global terrestrial carbon 
cycle model BETHY [Knorr, 1997] (see App. A). Thus, the variation of substrate availability 
with time t,f in (t), is parameterized as: 


/»«(*) ” 1 + Japp f NPp( l ) (2) 

r r max 

where NPP max denotes the annual maximum value of the NPP [gC*m' 2 *mo"']. In addition, the 
availability of potential substrate for methanogenesis decreases with increasing soil depth z, 
because the sources of fresh organic material are located either at the soil surface or in the 
upper soil layers. Since root exudates represent a major source of substrate, it is assumed that 
in vegetated soils the availability of substrate is constant throughout the rooting zone and then 
decreases exponentially with depth. The substrate availability f org (z) is then given by: 


forg(z) 


1 

(— \z - nroot\ )/( 10 cm) 

e 


nroot < z < ns 
nsoil < z < nroot 


(3) 


where nroot denotes the rooting depth, ns the soil surface and nsoil the soil depth, i.e. the lower 
boundary of the active layer. In unvegetated soils the vertical distribution of substrate f org (z) is 
assumed to decrease exponentially from the soil surface to the lowest soil layer: 

f org (z) = 0.857 . e CU-ns\)/(20cm) nsoil < z < ns (4) 
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The variation of the methane production rate with temperature is formulated in the following 
way. A Q 10 dependence is chosen using a Q l0 value of 6 lying within the range of observed 
QlO values ranging from 1.7 to 16 [Dunfield et al., 1993; Valentine et al., 1994; Westemiann, 
1993]. It is assumed that this temperature function includes both, the temperature dependence 
of the production of substrate for methanogenesis and that of methane production. The 
temperature function describes the response to the seasonal variation of the soil temperature 
T(t,z) at time t and depth z relative to the annual mean soil temperature T mean at the site. 
Consequently, the methane production rate R pro Jt,z) at time t and depth z is described as: 


T{t,z)-T m , on 

RproA 1 - Z) » *0 • /t,(2) ' /,„(<) ' f(T) ■ fiio 10 (5) 

R 0 is a constant rate factor expressed in units of pM/h (1M=1 mol/1). It is a measure for the 
.absolute* substrate availability and quality (since f in (t) and f org (z) describe only relative 
changes in the substrate availability with time and depth). Because in the model those 
quantities are not simulated explicitly, the parameter R 0 is a tuning parameter which has to be 
adjusted to each data set. As shown below, it only changes the amplitude but not the pattern of 
the modeled methane emissions. The function f(T) is a step function being 1, if T(t,z) is above 
0°C, and 0 else. This is equivalent to the assumption that there is no considerable methane 
production at sub-zero temperatures. There are a few studies that report winter methane 
emissions [Dise, 1992; Melloh and Crill , 1996] and at least one that reports methane 
production at temperatures down to -5°C [Clein and Schimel, 1995], but it does not seem that 
methane production at sub-zero temperatures is significant in general. In some of those studies 
it is also suggested that part of the methane emitted in the winter is methane which was 
produced already in the summer and was stored in the soil. Furthermore, there are several 
studies showing that methane emissions during the winter time are significantly lower than 
during the growing season [e.g. Whalen and Reeburgh, 1992; Shannon and White, 1994]. 

The chemical conditions in the soil can also affect methanogenesis [e.g. Conrad, 1989], For 
example pH conditions can be important. In addition, methanogenesis requires the absence of 
competing electron acceptors like sulfate (SO4 2 '), nitrate (N0 3 ‘) or ironlll (Fe 3+ ) and a redox 
potential below -200 mV. In addition, the availability of nutrients can enhance the productivity 
[Schlesinger, 1991]. Those effects are not taken into account in the model. The potential effect 
of a disadvantageous pH and the availability of nutrients on methane production rates is 
included in the parameter R 0 . Concerning the effect of competing electron acceptors and the 
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redox potential on methane production, it is assumed in the model, that the time scale, on 
which such effects are inhibiting methane production, are small compared to the time scale of 
the model, being one day. 

Methane oxidation rate R oxi( j 

The methane oxidation rate R ox jJt,z) at time t and depth z is assumed to follow Michaelis- 
Menten kinetics [Bender and Conrad , 1992], In addition, it depends on the soil temperature 
with observed Qiq values lying in the range between 1.4 and 2.1 [Dunfield et al., 1993; 
Knoblauch, 1994]. In the model, Q )0 for oxidation has been chosen to be 2. Thus, R ox jJt,z) is 
calculated from: 


R oxidi t > z) 


Yjnax ' ^CH4^ Z ) 
K m + z ) 


T(t,z)-T m 


Q 


10 


10 


(6) 


where K m and V max are the Michaelis-Menten coefficients. C CH4 (t,z) denotes the methane 
concentration [|iM] at time t and depth z, whereas T mean is the annual mean soil temperature 
[°C], The observed values for K m lie between 1 and 5 pM [Dunfield et al., 1993; Knoblauch, 
1994], while V max has been found to cover the range of about 5 to 50 pM/h [Dunfield et al., 
1993; Knoblauch, 1994; Krumholz et al., 1995; Moore and Dalva, 1997; Sundh et al., 1994; 
Watson et al., 1997], In the model K m and V max have been set to 5 pM and 20 pM/h, 
respectively. During the tests of the model against data sets from various wetlands it turned out 
that a V max value of 20 pM/h cannot be used at all sites. Therefore, at sites where a value of 20 
pM/h does not yield satisfactory results we choose another value for V max , the resulting values 
lying in the range from 3 to 45 pM/h. In the model, methane oxidation occurs only in the 
unsaturated soil layers. It has been observed, that in ecosystems with a water table being 
several decimeters above the soil surface there can be, due to turbulent diffusion, enough 
oxygen in the standing water column to make methane oxidation possible [Valentine, pers. 
comm.; Reeburgh, pers. comm.]. Since it is not known, if this effect is of general importance in 
different wetland types, it has not been incorporated into the model. 

Diffusion F di ff 

The diffusive flux F^Jt.z) is calculated using Fick‘s first law: 

F diff ( r > z) = • 0 ^^ c// 4 (r, z) (7) 
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where D CH4 (z) is the diffusion coefficient of methane at depth z and Ccttffhz) the methane 
concentration at time t and depth z. Since in the soil diffusion occurs only through the soil 
pores, the diffusion coefficient is obtained from: 

0Ctf4<z> = *>i - 0.66 /„„„ (8) 

which is the so-called Penman relation [ Hillel , 1982] using a tortuousity coefficient, which has 
been set to 0.66, suggesting that the distance covered by diffusion is about two-thirds of the 
length of the real average path. In the unsaturated soil layers, D i is the diffusion coefficient of 
methane in bulk air, which is 0.2 cm 2 /s [D'Ans and Lax, 1967] while in the water saturated soil 
layers D ( =10~ 4 *0.2 cm 2 /s [Scheffer and Schachtschabel, 1982], The factor f coarse denotes the 
relative volume of the coarse pores, i.e. we assume that diffusion proceeds mainly through the 
large pores. To solve the diffusion equation the Crank-Nicolson scheme [Press et al., 1992] is 
used. The methane concentration values are defined in the middle of each layer, whereas the 
diffusion coefficients are defined at the layer boundaries. The following boundary conditions 
are chosen: 



(/, z = nsoil) - 0 


(9) 


at the lower boundary nsoil and 

C CH4 (t,z = u+4cm) = C atm (10) 

at the upper boundary u, which is either the water table w(t) (if w(t)>ns) or the soil surface ns. 
C atm is the atmospheric methane concentration of 0.076 pM. At the water-air interface the 
methane concentrations in both phases are assumed to be in equilibrium. Hence their ratio is 
specified by the Bunsen solubility coefficient. 

Ebullition Q ebull 

As soon as the methane concentration in a layer exceeds a certain threshold concentration 
C thresh bubbles are formed. Since there are usually enough condensation nuclei in the soil and 
due to the occurrence of other gases such as nitrogen (N 2 ), C tres h usually is significantly lower 
than the saturation concentration (which is 1875 pM at 10 °C [Kaltofen et al., 1975]). In 
vegetated soils, for example, the concentration of N 2 was reported to be higher than in 
unvegetated soils [Chanton and Dacey, 1991]. Thus, according to Chanton and Dacey [1991], 
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Shannon et al. [1996] and Holzapfel-Pschorn et al. [1986] we assume that C tresh is lower at 
vegetated sites than at unvegetated sites. In the model it is assumed that bubble formation 
occurs at methane concentrations between 500 pM and 1000 pM at totally vegetated and 
unvegetated soils, respectively, being equivalent to a mixing ratio of (27-53)% of methane in 
the bubble (at 10°C). Thus, the threshold concentration for bubble formation is calculated 
from: 

C M = C min -{ l + %p) (11) 

where C min is the concentration at which bubble formation occurs, if the site is totally 
vegetated (500 pM) and P unveg is the percentage of unvegetated, bare soil. In order to 
determine the flux F ebu u(t ) of methane reaching the water table w(t) at time t in the form of gas 
bubbles, the rate Qetudhz) at which methane in the form of bubbles is removed from depth z is 
calculated: 


QebulM' k e ' f(C CH4 ) ■ (C CN4 (t, Z) C thresh ) (12) 

where k e is a rate constant of the unit 1/h and f(C CH4 ) is a step function taking the value 1, if 
the methane concentration C CH4 (t,z) is greater than C thresh , and 0 otherwise. Bubbles are 
assumed to reach the water table within one model timestep being 1 hour. This assumption is 
supported by the facts that wetland soils are generally very porous, the relative pore volume 
being often greater than 90% [Scheffer and Schachtschabel, 1982], and the finding that the 
velocity of bubbles ascending in pure water lies in the order of 1-10 cm/s [e.g. Shafer and Zare, 
1991]. Thus, the ebullitive flux F ebu n(t) is obtained by integrating Q ebu u(t,z) over the whole 
water saturated zone: 


F ebM (i) = r (0 ,e rf „„(<.z)* 03) 

J nsoil 

where nsoil is the lower boundary of the active layer and w(t) the water table. If the water table 
is at or above the soil surface, F ebull (t) contributes directly to the methane flux into the 
atmosphere. Otherwise, the amount of methane rising to the water table in the form of bubbles 
is added to the methane concentration in the lowest unsaturated soil layer. 

Plant-mediated transport Q p i ant 

The occurrence of vascular plants provides an effective mechanism by which methane can be 
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transported up to the atmosphere [e.g. Chanton et al., 1992a; Schimel, 1995; Shannon et al., 
1996]. Moreover, oxygen can be conducted down to the roots this way establishing a small 
aerobic region in the rhizosphere [Gerard and Chanton, 1993; Holzapfel-Pschom et al., 1986; 
Schipper and Reddy, 1996; Schiitz et al., 1989]. The results from several studies of plant- 
mediated transport suggest that the main emission pathway is by molecular diffusion or 
effusion through the plant stems [ Chanton et al., 1992a; Chanton et al., 1992b; Happell et al., 
1993; Nouchi and Mariko, 1993; Shannon et al., 1996]. Consequently, the flux of methane 
through plants is assumed to be proportional to the concentration gradient between the soil and 
the atmosphere and - since the atmospheric methane concentration is very small compared to 
the methane concentration in the soil - proportional to the methane concentration in the soil. 
This assumption corresponds well with observations by Nouchi and Mariko [1993] who found 
a linear relationship between methane emission rates through plants and pore water methane 
concentrations. On the other hand, methane entering the roots of plants has to pass through the 
small oxic zone around the root tips and a certain fraction P ox of methane is oxidized. There 
are a few studies examinig rhizospheric oxidation: Schipper and Reddy [1996] determined the 
fraction of methane consumed by rhizospheric oxidation using two different methods and 
found values of 65±24 % and 79±20%. Schiitz et al. [1989] observed that about 90% of the 
methane produced was oxidized in case plant-mediated transport was the main (97%) emission 
pathway, while Gerard and Chanton [1993] obtained values for rhizospheric oxidation 
between 39% and 98%. Taking the high variability of observed values we set P ox to 50%. 
Hence, the rate Q p i ant (t,z) at which methane is removed by plants from depth z at time t is 
calculated from: 

Qplaruit’Z) = - k p r veg- f rool^) ' f grow(0 ' C CH4 (t, Z ) ■ (l - P ox ) ( 14 > 

where k p is a rate constant of the unit 0.01/h and T veg a factor describing the , quality* of plant- 
mediated transport at a site, depending on the density of plant stands and the plant types and 
being in the range from 0 to 15. At sites where the predominant plant types and their capability 
of conducting gas are known, T veg is derived from that knowledge. Otherwise, we consider 
shrubs not contributing to plant-mediated transport and trees being poor, grasses and sedges 
being good gas transporters. The function f roo /z) represents the vertical distribution of roots in 
the soil. It is assumed to decrease with depth and for the sake of simplicity we choose a linear 
relationship between root biomass and soil depth: 


\^-nroot\ + 1 
k root ns-nroot 


(15) 


10 



where nroot denotes the rooting depth and ns the soil surface. The ability of plants to conduct 
gas is considered to vary with the growing state of the plants. This assumption agrees well with 
observations by Schiitz et al. [1989], who found that the fraction of methane transported 
through rice plants increases with growing maturity of the plants. In the model we use the leaf 
area index (LAI) as a measure for the growing state f grow (t) of the plants (see App. B). Finally, 
the methane flux due to plant-mediated transport F p i an j(t,z) is calculated from: 

F plml U) = J" Q pl „„«)* (16) 

1 • n rnnt 1 


where ns and nroot are the soil surface and the rooting depth, respectively. 

Total methane emission F tot 

The total methane emission F tot (t) is calculated by adding all the fluxes from the different 
transport mechanisms. The diffusive flux F di ^t,z=u) at the soil/water-atmosphere boundary u, 
the flux due to plant-mediated transport F plan Jt,z ) and the ebullitive flux F ebull {t,z). The latter 
contributes only to the total flux F tot (t), if the water table is at or above the soil surface. Hence 
F tot (t) is obtained from: 


F , o ,(0 = F diff( t ’ z = U ) + F ebuii( t )+F planl (t) (17) 


Results and Discussion 

The methane model is tested against observational data from 15 microsites within 6 wetlands 
located in North America, Europe and Central America. The observation periods range from 1 
season up to 3 years. Tab. 1 gives a summary of the data sets used. At each site methane 
emissions were measured about 2-4 times per month and at two sites (1, 2) also methane 
concentration profiles in the soil were obtained 1-2 times per month. At all sites except site 5 
the position of the water table and the soil temperature at different soil depths were observed at 
least at the same frequency as the methane emissions. Since at site 5 the forcing data were not 
observed a hydrologic model [Walter, 1998] in combination with the ECHAM4 model 
[Roeckner et al., 1996] is used to run the methane model there. The daily forcing required by 
the methane model is obtained by linear interpolation of the observed input data. The NPP is 
obtained from the BETHY model [Knorr, 1997] for that 1°* 1° grid box in which the wetland is 
located. The parameters soil depth nsoil, rooting depth nroot and the quality of plant-mediated 


11 



transport T veg are chosen based on information about the specific wetland site given by the 
investigators. The parameter V max is set to 20 pM/h in most cases, but is adjusted to the data set 
in situations where a value of 20 pM/h does not yield optimal results. The values of V max used 
in this study range from 3 to 45 pM/h. The parameter Rq is adjusted to the data set in each case. 
The values for R 0 obtained this way lie between 0.3 and 0.6 at the high latitude sites and is 2.8 
at the tropical site. In this article results of 5 tests from one station at every site 1-5 are shown. 
The results from all 15 microsites are compiled in detail in Walter [1998]. Tab. 2 lists the 
parameters used at the 5 test sites. 


Site 1: Michigan 

Fig. 2 shows the results of a test of the model at site 1 using data by Shannon and White 
[1994], The forcing data are plotted in the lower part, (c) shows the observed position of the 
water table relative to the soil surface and (d) the observed soil temperatures at different soil 
depths. In (a) the comparison between simulated and observed (the average of 3 chamber 
measurements ± 1 standard deviation error bars) methane emissions is shown. When the water 
table is above the soil surface the patterns of observed as well as simulated methane emissions 
are dominated by changes in the soil temperature. This applies to both day-to-day and 
interannual variations of methane emissions. For example, methane emissions are higher in 
1991 than in the 2 following years, because of slightly higher soil temperatures in 1991. As 
soon as the water table falls below the soil surface this connection between methane emission 
and soil temperature is no longer valid. In this situation methane diffusing through the soil 
pores is partly oxidized in the unsaturated soil layers. In addition, rising bubbles cannot reach 
the soil surface anymore and only the deeper roots extend into the water saturated soil zone 
where high methane concentrations prevail. The consequence of a long period with the water 
table being below the soil surface can be seen in 1991, when methane emissions drop to 
relatively small values of about 200 mg^m^^d' 1 already in August even though the soil 
temperature is still high at that time. The reason for this is that the water table stays at depths of 
about 15 cm below ground since July 1991. In 1992 and 1993 this drop of methane emissions 
occurs much later in the year, attributable to wetter conditions. In (b) the simulated relative 
contributions of the three different transport mechanisms diffusion, ebullition and plant- 
mediated transport are depicted. At sites covered by vascular plants, as it is the case at site 1, 
plant-mediated transport plays an important role during the growing season. Shannon et al. 
[1996] observed in plant enclosure experiments that 64-90% of the net methane flux was 
emitted through Scheuchzeria palustris. The model results are in good agreement with this 
finding (Fig. 2 (b)). Ebullition occurs only at times when the water table is above the soil 
surface. Diffusion plays only a role if the water table is below the soil surface. For example. 


Figure 2 
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immediately after the fall of the water table below the soil surface in June 1991 diffusion 
increases considerably, because the diffusion coefficient increases by a factor of 10 4 as soon as 
the soil becomes unsaturated and methane concentrations are still high in the uppermost layers. 
After a period of about 15 days, emission by diffusion declines again due to oxidation in the 
unsaturated zone. Moore and Roulet [1993] observed a similar pattern: increased fluxes with 
falling water table to 20 cm depth within 10 days, followed by decreased fluxes as the water 
table continues to fall further. 


Fig. 3 shows a comparison between simulated and observed methane concentration profiles for 
the period between January 1992 and January 1993. The observations as well as the model 
results show the same seasonal pattern: higher methane concentrations in the winter and spring 
and decreased methane concentrations in the summer and autumn. The reason for this is that 
during the growing season a large fraction of methane is removed from the soil by plants. 
However, in the model removal of methane from the soil by plants seems to start too early 
(May 16 and Jun 4). This probably happens, because in the model the time plants need to 
develop is underestimated. At the beginning of the winter the concentrations increase faster in 
the data than in the model (Nov 14, Dec 13 and Jan 9), again because the simulation of the 
seasonal variation of the growing stage of the plants seems to fail in late autumn. However, 
these discrepancies occur only at times when the production rate is relatively low, and thus 
affect simulated methane emissions only slightly. 


Figure 3 


Site 2: Minnesota 

Fig. 4 shows the results of a model-data comparison performed at site 2 with data from Dise 
[1993]. As above, in the lower part the observed forcing data which are the position of the 
water table (b) and the soil temperature (c) are plotted. In (a) the comparison between observed 
and simulated methane emissions is shown. Here, the agreement seems to be less good than at 
site 1 which can be partly explained by the fact that at site 1 the model results are compared to 
an average of 3 measurements, which clearly differ from each other at some times (see 
±1 standard deviation error bars in Fig. 2 (a)), whereas at site 2 only 1 chamber was used. In 
addition to that, other factors affecting methane production and emission such as the chemical 
conditions in the soil, which are not considered in the model (see above), could be important at 
some times. Generally, simulated methane emissions respond in a similar way to changes in 
the soil temperature and the position of the water table as at site 1. 


Figure 4 


Site 3: Finland 

In Fig. 5 the results of a comparison between modeled and observed methane emissions from 2 
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microsites (flark and low hummock) of a Finnish wetland [Saamio etal., 1997] are shown. We 
focus on intrasite variations within the same wetland, here on differences in the position of the 

water table. The same values for R 0 , V max , nsoil, nroot and T veg are used at the two microsites. 

In addition, the soil temperature as shown in Fig. 5 (c) is the same at both microsites. Hence Figure 5 
the differences in simulated methane emissions between these two microsites are attributed to 
differences in the positions of the water table only. The low hummock is elevated by about 15 
cm relative to the wetland surface. Fig. 5 (b) shows the positions of the water table at the low 
hummock and the flark. The modeled and observed methane emissions from these 2 microsites 
are plotted in Fig. 5 (a) (the grey areas are averages of the fluxes from all the collars of the two 
microsites +1 standard deviation). The methane emissions from the flark are considerably 
higher than those from the low hummock, because there is less oxidation due to the higher 
water table level at the flark site. This effect can be seen both in the observations and the model 
results. At the flark simulated methane emissions occur slightly too early (see Fig. 5 (a)). The 
modeled fluxes follow directly the changes in the soil temperature (Fig. 5 (c)) whereas in the 
observations there is a time lag between the soil temperature and methane emissions. The 
reason for this time lag could be that before methane production can start, appropriate soil 
conditions for methanogenesis must be established and suitable substrate for methanogenesis 
must be produced. This is not taken into account in the model. Another reason could be that at 
these microsites the deeper soil layers contribute more to methane production than the upper 
ones, possibly due to better substrate availability and/or quality. This would cause modeled 
methane emissions to follow the soil temperature of deeper soil layers which lags behind the 
temperature of the upper soil layers. However, this effect is only seen at that two microsites of 
the Finnish wetland. 

Site 4: Alaska 

A further test is performed with data from an Alaskan Arctic tundra site which is underlain by 
permafrost [Whalen and Reeburgh , 1992]. In Fig. 6 (b) and (c) the observed forcing data water Figure 6 
table, thaw depth and soil temperature are plotted. Since the three chamber measurements at 
the black hole microsites (called BH1-3) differ much, we plot the obtained fluxes from each 
chamber as well as the average of these measurements to compare them with the simulated 
fluxes (Fig. 6 (a)). The vegetation consists of Sphagnum mosses and there are no vascular 
plants which means that plant-mediated transport does not occur and hence T veg = 0 (see Tab. 2). 

During the summer the soil starts to thaw and the maximum thaw depth of about 60 cm below 
ground is attained in autumn. An increasing thaw depth means that the production zone is 
being enlarged with time. Therefore, simulated methane emissions are of the same magnitude 
in September as in August even though both the soil temperature and the water table levels are 
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lower in September. Here V max , which describes the maximum methane oxidation rate in the 
unsaturated soil zone, has been set to only 3)tM/h as opposed to values of 20pM/h at most 
other sites (see Tab. 2). Because of that low V max value, the effect of a declining water table on 
methane emission is smaller than at the Finnish wetland site discussed in the previous 
paragraph. One could speculate what the reason for this low V max value might be. For example, 
the oxidation potential of the methanotrophic bacteria could be lower at that site or the soil 
moisture content in the soil layers above the water table could be higher (the soil moisture in 
the unsaturated zone is not taken into account in the model). However, the V max value chosen 
here is still close to observed values for V max lying in the range between 5 and 50 |iM/h 
[Dunfield et al., 1993; Knoblauch, 1994; Krumholz et al., 1995; Moore and Dalva, 1997; 
Sundh et al., 1994; Watson et al., 1997]. 


Site 5: Panama 

As a last example a test against data from a tropical wetland in Central Panama [Keller, 1990] 
is discussed. At this site methane emissions were measured, but neither the soil temperature 
nor the water table was observed. This site is characterized by a 4-month dry season between 
February and May. According to Keller [pers. comm.] the water table rarely exceeded +30 cm 
during the wet season and was below 50 cm below ground during the dry season, and the soil 
temperature was nearly constant over the whole year. In order to obtain input data for the 
methane model for this test the output of a climatological run of the ECHAM4 model 
[Roeckner et al., 1996] in T 106 (1.1 °* 1.1°) is used. The soil temperature at different soil 
depths is calculated by the ECHAM4 model and the water table levels are calculated by a 
hydrologic model [Walter, 1998] driven by climate data from the same run of the ECHAM4 
model. Since the necessary model parameters (see Tab. 2) are not known for this site the values 
from global data sets of these parameters which have been established for a global application 
of the methane model [Walter, 1998] are used. In Fig. 7 (b) the water table calculated by the 
hydrologic model forced by the ECHAM4 model is plotted. The soil temperature is almost 
constant throughout the year (not shown). This is consistent with the information by Keller 
[pers. comm.] (see above). Fig. 7 (a) shows the comparison between the simulated and the 
observed methane emissions. The latter are averages of 7 stations plotted with ±1 standard 
deviation error bars. The most striking point is that the modeled as well as the observed 
methane emissions are zero or even negative during the dry season, attributable to the decline 
of the water table to depths of about 80 to 100 cm below ground, causing the development of a 
large oxic zone where methane is consumed. This data set is the only tropical data set where 
methane emissions were measured over a period of at least 1 season that has been available to 
us until now. Since it is possible that in tropical wetlands processes are dominant which are not 


Figure 7 
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yet included in the model, additional data sets from tropical wetlands are needed, in order to 
further improve and test the model. This is particulary important, because tropical wetlands are 
a very larger methane source [e.g. Hein et al., 1997]. 


Relationship between climate input and methane emissions 

The time series used above to test the model show that the pattern of methane emissions at 
most high latitude sites (1-3) is mainly controlled by soil temperature as long as the water table 
is above the soil surface. If the water table is below the soil surface and does not vary much, 
like at site 3 (Fig. 5, hummock), the pattern of methane emissions is again mainly driven by 
changes in soil temperature, but the amplitude is influenced by the water table depth. However, 
as soon as the water table fluctuates around the soil surface, as for example at site 1, this 
connection between soil temperature and methane flux does no longer apply. As it is shown in 
Fig. 2 and discussed above methane emissions increase in the first about 15 days after a fall of 
the water table below the soil surface. This behavior is also reflected in Fig. 8 which shows the 
connection between simulated methane emissions and soil temperature (a) and water table (b) 
at site 1. Fig. 8 (a) shows that there is a positive connection between soil temperature and 
simulated methane emissions. However, in the temperature range between 15 and 20°C, for 
example, simulated methane emissions range between 50 and 900 mg*m" *d' . Fig. 8 (b) 
shows that there is no connection between water table and simulated methane emissions at site 
1. The highest methane emissions occur, when the water table is 5cm above or 5cm below the 
soil surface. When the water table is 13 or 17cm below ground, methane fluxes can be as high 
as 450 and 250 mg*m' 2 *d 1 , respectively. This agrees well with the results of Shannon and 
White [1994] who calculated the correlation between soil temperature and water table and 
observed methane fluxes at that site. At the Arctic tundra site (site 4) the connection between 
the environmental parameters soil temperature, water table and thaw depth and simulated 
methane emissions looks different (Fig. 9). There is no obvious connection between soil 
temperature and simulated methane emissions (a) and the same applies to the connection 
between water table depth and simulated methane emissions (b). However, simulated methane 
fluxes are higher at lower thaw depths (b), i.e. when a larger part of the soil is thawed. Whalen 
and Reeburgh [1992] found that methane emissions from all stations investigated at that Arctic 
tundra site correlated best with ‘centimeter-degrees’, being defined as the absolute value of the 
product of thaw depth and mean soil temperature to permafrost. At the tropical test site (site 5) 
the temperature does not vary much and from Fig. 7 it is clear that their is a connection 
between the position of the water table and simulated methane emissions. These results show 
that in most cases simple relationships between environmental parameters and methane fluxes 
do not even apply at one particular site. They differ from site to site and hence a relationship 


Figure 8 
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between climate input and methane emissions found at one site cannot be applied to other sites 
and other conditions. This means that if one wants to model methane emissions from different 
wetlands and under varying environmental conditions, a more process-based approach is 
needed. 


Sensitivity Tests 


Sensitivity to model parameters 

The sensitivity of the model to the choice of the parameters listed in Tab.2 is tested at site 1 . 
The only other site-specific model parameter not listed in Tab. 2 is the relative volume of 
coarse pores f coarse , which has been set to 0.5 at all sites except site 5, where it has been 
derived from a global data set used for the global model run as discussed above and in Walter 
[1998]. The model is only sensitive the choice of f coarse at sites where diffusion is the dominant 
transport mechanism which is only the case at site 4 which is unvegetated and where the water 
table is below the soil surface most of the time. A sensitivity test carried out at site 4 shows that 
a change in f coarse of +0.1/-0. 1 leads to a 12.5% increase/decrease in simulated methane 
emissions at that site (not shown). Differences in the soil depth nsoil only rarely affect 
simulated methane emissions (provided nsoil<nroot; not shown), because below the rooting 
depth the availability of substrate for methane production decreases exponentially (see Eq. 
(3)). In Fig. 10 the results of the sensitivity tests with R 0 , T veg and moot are shown. The 
results of the test with R 0 (see Eq. (5)), which is adjusted to each data set, show that R 0 
influences the amplitude of simulated methane emissions, but not the pattern (Fig. 10 (a)). The 
choice of the parameter V max changes the model results only in situations when the water table 
is below the soil surface (Fig. 10 (b)). In this situation a lower V max leads to higher methane 
emissions, because of lower methane oxidation rates (see Eq. (6)). The parameter T veg 
influences the fraction of methane which is transported through plants, but also the amount of 
methane that is oxidized in the rhizosphere (see Eq. (14)). Its influence on simulated methane 
emissions depends on the water table level. If the water table is above the soil surface a higher 
T veg leads to lower methane emissions (Fig. 10 (c)), because in this case less methane is 
transported by ebullition and more through plants. Therefore, more methane is oxidized in the 
rhizosphere. This agrees well with observations by Holzapfel-Pschom etal. [1986] who found 
that there was more ebullition at unvegetated sites than at vegetated ones and that a larger 
fraction of produced methane was emitted from unvegetated sites, if the water table was above 
the soil surface. If the water table is below the soil surface, diffusion is high in the first month 
after the fall of the water table below the soil surface (see Fig. 2 (b)). Since a lower T veg leads 
to higher diffusion rates, because less methane is removed from the soil in this case (see Fig. 

1 1), a lower T veg leads to higher methane emissions in the first month after the water table fell 


Figure 10 
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below the soil surface. In situations where the water table is below the soil surface for a longer 
period of time a higher T veg leads to higher methane emissions, because then a huge fraction of 
methane diffusing through the oxic soil layer is oxidized and methane transported through 
plants bypasses this oxic layer [Walter, 1998]. In Fig. 1 1 the effect of T veg on the simulated soil 
methane concentration is shown. It is clear, that higher T veg values lead to lower soil methane 
concentrations, because more methane is removed from the soil by plants. This agrees well 
with observations that the soil methane concentration is higher at unvegetated sites than at 
vegetated sites [Chanton and Dacey, 1991]. The rooting depth nroot is the soil depth down to 
which plants can extract methane from the soil (see Eq. (16)). Furthermore, nroot affects the 
vertical distribution of substrate in the soil (see Eq. (3)). Throughout the root zone substrate 
availability is constant, then it decreases exponentially. In the model the combined effect of an 
enlarged production zone and more plant-mediated transport and hence more rhizospheric 
oxidation are almost balanced, if the water table is above the soil surface (Fig. 10 (d)). 
However, if the water table falls below the soil surface transforming the soil above the water 
table from a production into a consumption zone, the enlargement of the zone with high 
productivity (higher nroot ) shows a pronounced effect on simulated methane emissions. 


Figure 11 


Sensitivity to changes in the input data (climate forcing) 

The effect of changes in the soil temperature and water table on simulated methane emissions 
are shown in Fig. 12. The effect of the water table being raised or lowered by 10cm (Fig. 12 
(b)) on simulated methane emissions is plotted in Fig. 12 (a). As long as the water table is 
above the soil surface the height of the water table does not change simulated methane 
emissions. If the water table is below the soil surface methane emissions decrease with 
decreasing water table. Therefore, the effect of a climatic change leading to a shift in the 
position of the water table at a given site depends on the conditions prevailing at that site. If the 
site is relatively dry like in 1991 a shift in the position of the water table is very pronounced 
(see Tab. 3), whereas if the site is relatively wet as in 1993 the same shift in the water table 
changes methane emissions much less. An uniform change in the soil temperature by ±1°C 
leads to a regular shift in the amplitude of simulated methane emissions which does not depend 
on the position of the water table (Fig. 12 (c)). An uniform increase in the soil temperature of 
1°C leads to an increase in simulated methane emissions of about 20.5% at that site in all three 
years (see Tab. 3), whereas an uniform decrease in the soil temperature of 1°C reduces 
methane emissions by about 17%. These sensitivity studies demonstrate that possible climatic 
change in the future could considerably alter methane emissions from natural wetlands. 


Figure 12 


18 




Summary and Conclusion 

Methane emissions from natural wetlands, which constitute the biggest natural methane source 
at present, are highly climate-sensitive. In this study we presented a model to simulate the 
processes leading to methane emission within a 1 -dimensional soil column as a function of 
climate. The model differs from all other models in the literature in the way that transport is 
parameterized. Two opposing processes operate in soil, namely methane production in the 
anoxic soil zone and methane consumption in the oxic soil zone. Therefore transport, which 
occurs by diffusion, ebullition and through plants, plays an important role in determining the 
fraction of produced methane that is emitted into the atmosphere. Moreover, the different 
transport mechanisms influence the velocity at which methane is transported to the 
atmosphere. For climate forcing variables soil temperature, water table and thaw depth (which 
had been observed at 4 sites and were simulated at one site) as well as model derived NPP were 
used. Site-specific model parameters are soil depth nsoil, rooting depth nroot, the relative 
volume of the coarse pores f coarse , the quality of plant-mediated transport T veg , and in some 
cases the maximum methane oxidation rate V max . The chemical properties (including substrate 
availability and quality for methanogenesis) of the soil and their influence on methane 
production are not modelled explicitly, but are included in the parameter R 0 which has been 
adjusted to each data set. In order to apply the model to other sites without tuning R 0 or to 
larger (regional or global) scales, a method must be developed to derive R 0 from 
biogeochemical, biogeographical and climatic variables. In the future a process-based model to 
simulate the production of substrate for methanogenesis as a function of the soil conditions and 
the climate could also be used to replace the tuning parameter R 0 . The occurrence of turbulent 
diffusion in the standing water and its effect on the transport of methane to the atmosphere and 
re-oxidation of methane in the water column could be important at tropical sites in particular. 
However, it has not yet been included in the model, as the data sets used in this study showed 
no evidence that this process is of general importance. 

In this article results of tests of the model against observational data from 5 different wetland 
sites located in North America, northern Europe and Central America, representing a large 
variety of environmental conditions, are presented. The observational periods range from one 
season to several years. One site is underlain by permafrost and undergoes a seasonal freeze 
thaw cycle, while one site is located in the tropics where there is a dry/wet-seasonal cycle. At 
one site simulated methane concentration profiles in the soil were compared with observed 
vertical concentration profiles which reflect the vertical distribution of methane production, 
methane oxidation and removal of methane by the different transport mechanisms. Hence, 
methane concentration profiles constitute an additional constraint to the model. We 
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investigated the climate-induced variations of simulated methane emissions at the test sites and 
conclude that this model can be applied to different wetlands under various conditions and can 
also be applied to the global scale. The relationship between the input data (soil temperature, 
water table and thaw depth at the permafrost site) and simulated methane emissions was 
examined. Simple relationships between the input data and simulated methane emissions do 
not, in general, apply at single sites and they differ from site to site. This means that a more 
process-based approach is needed in order to simulate climate-induced variations of methane 
emissions from different wetland sites. Sensitivity tests of the model to the choice of model 
parameters were performed. The tuning parameter Rq (in the methane production rate 
equation) changes the amplitude, but not the pattern of simulated methane emissions. The 
choice of the parameter V max (in the methane oxidation rate equation) alters simulated methane 
emissions in situations when the water table is below the soil surface. The response to different 
choices of the parameter T veg (quality of plant-mediated transport) depends on the fact, if the 
water table is above or below the soil surface and in the latter case for how long. In addition, 
the sensitivity of the model to uniform changes in the input data (soil temperature, water table) 
was tested. The response of a uniform change in the soil temperature of +1°C or -1°C is in the 
order of a 21% increase/ 17% decrease in methane emissions, respectively, no matter if the soil 
if relatively dry or wet. However, the response of the model to uniform changes in the water 
table of ±10cm depends highly on the prevailing hydrologic conditions in the soil. In a 
relatively dry year an increase in the water table of 10cm can enhance simulated methane 
emissions by as much as 25%, whereas in a wet year there can be no change at all. Similarly, a 
decrease in the water table of 10cm affects simulated methane emissions much more in a 
relatively dry year (in this study 37% less methane is emitted) than in wetter years. These 
numbers suggest that possible climate change in the future could considerably affect methane 
emissions from natural wetlands. However, in order to assess the change in methane emissions 
due to future possible climate change a more realistic scenario simulation should be carried 
out. It should use simulated future climate as input data and include possible substrate 
limitation, changes in wetland areas and melting of permafrost. The latter is included in the 
model via the soil temperature, since the soil layers where the soil temperature is below 0°C 
are considered to be frozen. 

Further studies and data sets are necessary to further improve the model and test it more 
thoroughly. For example, only one data set from tropical wetlands was available to us. Since 
tropical wetlands are a very important methane source, we propose that in the future emphasis 
should be placed on long-term studies of methane emissions from tropical wetlands. In 
addition, it would be advantageous to test the model with a data set consisting not only of time 
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series of the input and output data of the model (as it was the case in this study), but also of 
quantities calculated in the model such as methane production rate, methane oxidation rate and 
the fraction of methane transported by the different transport mechanisms. No such data set is 
currently available and the development of such a data set is therefore a priority to enable a 
more thorough testing of this model. 
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Appendix A: Calculation of f NPP (t)\ 


The function f^ PP (t) is derived from the variation of the NPP with time, NPP(t), as calculated 
by the BETHY (Biosphere-Energy Transfer and Hydrology) model [Knorr, 1997], The 
BETHY model is a process-based model describing the water balance on vegetated surfaces 
and bare soils and the C0 2 balance in vegetation and soils. It uses remote sensing data and 
calculates the Net Primary Productivity (NPP) on a 0.5°*0.5° grid with monthly time steps. 
NPP(t) is used as a measure for the variation of substrate availability with time. Since part of 
the substrate is degraded from organic matter incorporated into the soil from dying plants in 
autumn, we presume that in regions with a change between growing and non-growing season, 
i.e. where the growing season lasts for 3-9 months, the substrate availability increases again in 
autumn. We define the duration of the growing season as the time span when the soil 
temperature at 50cm depth below ground is above 5°C. 


In regions where the growing season is shorter than 3 months or longer than 9 months f^ PP ( t) 
is: 

f NPP (t) = NPP(t) at any time (18) 

whereas in regions with a growing season lasting 3 to 9 months f^ PP (t) is obtained from: 


NPP(t) 


f NPpW ~ 


NPP mn -NPP. t 

fsjpp _i max last ^ 

last 0 5 -/ 

nongrow 


~ NPP r , r „ - 

V ”"“l 0.5 • , ntmgrow 


growing season 

1. half non-growing season 

(19) 

2. half non-growing season 


Here, t nongrow is the duration of the non-growing season in days, and NPP tas[ and NPPjj rst 
denote the NPP(t) values of the last and first day of the growing season, respectively. 


22 



Appendix B: Calculation of f grow (t): 


The function f grow (t) describing the growing state of the plants is assumed to be a function of 
the LAI and is calculated according to Dickinson et al. [1993]. We define a mean LAI 
depending on the daily mean temperature at 50cm depth below ground T 5 o= T(t,z= 120cm) and 
hence calculate f grow (t) from: 


^50 < T grow 

T grow^T 5 10 <C, (20) 

T mat < Tso 

where A min and A have been chosen to be 0 and 4, respectively, and A max =A min +A. T grow is the 
temperature at which plants start to grow. In regions where the annual mean soil temperature is 
below 5°C, T grow has been set to 2°C, elsewhere to 7°C. This distinction has been made, 
because in cold regions plants start to grow at lower temperatures than in temperate and warm 
regions. The same is valid for the temperature at which they reach maturity, T mat , which has 
been set to T mat =T grow + 10°C. 


f grow(^) 


A min + A ' 1 1 “ 


\ A max 


T mat ~ ^50 

T~ -T 

mat growJ 
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Figure Captions: 


Fig- 1 : Schematic representation of the model structure: the 1 -dimensional soil column is divided 
into 1cm thick parallel layers. The forcing consists of daily records of the water table, the soil 
temperature and the NPP. Methane production occurs in the soil zone between the soil depth nsoil 
and the water table w, which can be either below or above the soil surface ns. The methane pro- 
duction rate R pro Jt,z) is a function of the soil temperature T soil (t,z) and the NPP, which is taken as 
a measure for substrate availability. Methane oxidation is confined to the soil layers between the 
water table and the soil surface. The methane oxidation rate R ox iJ.t,z) follows Michaelis-Menten 
kinetics and is a function of the soil temperature T soi [(t,z). Transport proceeds by (1) molecular 
diffusion through the soil pores, (2) ebullition, which is the formation of gas bubbles in the water 
saturated layers and their ascent to the water table and (3) plant-mediated transport from layers 
above the rooting depth nroot to the atmosphere. The model calculates methane fluxes to the 
atmosphere and methane concentration profiles in the soil on a daily basis. 

Fig. 2: Site 1: (a): comparison between modeled (thick line) and measured (dots with +/- 1SD 
error bars) methane emissions, (b): modeled contributions of the three transport mechanisms: dif- 
fusion (black), ebullition (light grey), and plant-mediated transport (dark grey), (c): forcing: 
observed position of the water table, (d): forcing: observed soil temperature at different soil 
depths, (observational data from Shannon and White [1994]) 

Fig- 3: Site 1: comparison between modeled (thick line) and observed (squares) methane concen- 
trations in the soil, (observational data from Shannon and White [1994]) 

Fig. 4: Site 2: (a): comparison between modeled (thick line) and measured (dots) methane emis- 
sions. (b): forcing: observed position of the water table, (c): forcing: observed soil temperature at 
different soil depths, (observational data from Dise [1993]) 

Fig. 5: Site 3: comparison of methane emissions from two microsites, (a): modeled fluxes from 
the flark (solid line) and the low hummock (dot-dashed line); the shaded areas are the areas 
between the average of all measurements from the respective microsite +/- 1 SD: flark (dark grey) 
and low hummock (light grey), (b): forcing: observed positions of the water table at the flark 
(solid line) and the low hummock (dot-dashed line) relative to the soil surface (dashed line), (c): 
forcing: observed soil temperature at different soil depths (the same at both microsites), (observa- 
tional data from Saamio et al. [1997]) 



Fig. 6: Site 4: comparison between modeled (thick line) and measured (the triangles and the 
square are observations from three stations within the , Black Holes 1 site, the dashed line marks 
the average of those three measurements) methane emissions, (b): forcing: observed thaw depth 
(soild line) and water table (dot-dashed line), (c): forcing: observed soil temperature at different 
soil depths, (observational data from Whalen and Reeburgh [1992]) 

Fig. 7: Site 5: (a): comparison between modeled (thick line) and measured (dots with +/- 1SD 
error bars) methane emissions, (b): forcing: position of the water table relative to the soil surface 
as calculated with the hydrologic model forced by the ECHAM4 model (see text), (observational 
data from Keller [1990]) 

Fig. 8: Site 1: relationship between input data and simulated CH 4 emissions: (a): connection 
between the soil temperature T so jj at the soil surface and simulated CH 4 emissions, (b): connec- 
tion between the position of the water table and simulated CH 4 emissions. 

Fig. 9: Site 4: relationship between input data and simulated CH 4 emissions: (a): connection 
between the soil temperature T so j| at the soil surface and simulated CH 4 emissions, (b): connec- 
tion between the position of the water table (grey circles) and the thaw depth (black triangles) and 
simulated CH 4 emissions. 

Fig. 10: Sensitivity to R 0 , V max ,T veg and nroot tested at sitel: (a): simulated CH 4 emissions using 
three different values for R 0 , namely 0.6 (black, soild), 0.3 (grey) and 0.9 (black, dashed), (b): 
simulated CH 4 emissions using three different values for namely 45pM/h (black, soild), 
20pM/h (grey) and 3pM/h (black, dashed), (c): simulated CH 4 emissions using three different 
values for T veg , namely 15 (black, solid), 5 (grey) and 1 (black, dashed), (d): simulated CH 4 emis- 
sions using three different values for nroot, namely 50cm (black, solid), 30cm (grey) and 70cm 
(black, dashed). 

Fig. 11: Sensitivity to 7 tested at site 1: simulated CH 4 concentration profiles in the soil using 
three different values for T veg , namely T veg = 15 (black, solid), T veg =5 (grey) and T veg = 1 (black, 
dashed) and comparison with observations (squares). 

Fig. 12: Sensitivity to forcing data (water table (a, b) and soil temperature (c)) tested at site 1 : (a): 
simulated CH 4 emissions using the original water table (black, solid), a by 10cm lowered one 
(grey) and a by 10cm raised one (black, dashed), (b): forcing: different water table levels relative 



to the soil surface (grey, dashed) used to test the sensitivity of the model to the water table: origi- 
nal water table (black, solid), by 10cm lowered water table (grey, soild), by 10cm raised water 
table (black, dashed), (c): simulated CH 4 emissions using the original soil temperature (black, 
solid), a by 1°C lowered one (grey) and a by 1° raised one (black, dashed). 
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Table Captions: 


Tab. 1: Locations, numbers of investigated stations and references of the data sets used to test 
the model. 

Tab. 2: Model parameters used at the 5 test sites presented in this article: R 0 (tuning 
parameter), V max (maximum methane oxidation rate), nsoil (soil depth), nroot (rooting depth) 
and T veg (quality of plant-mediated transport) (see model description). 

Tab. 3: Results of the sensitivity tests to the input data. Row 1-2: in runs T+1°C/T-1°C the soil 
temperatur was uniformly increased/decreased by 1°C in all soil layers throughout the whole 
simulation period. Row 3-4: in runs wt+10cm/wt-10cm the position of the water table was 
uniformly increased/decreased by 10cm throughout the whole simulation period (see Fig. 12). 
Row 5: simulated annual methane emissions for the years 1991-93. 



Table 1: Data sets used to test the model 


No. 

Site 

Location 

Stations 

Reference 

1 

Michigan 

42°N, 84°W 

1 

Shannon and White [1994] 

2 

Minnesota 

47°N, 93°W 

5 

Dise [1993] 

3 

Finland 

63°N, 31°E 

4 

Saamio et al. [1997] 

m 

Alaska 

65°N, 148°W 


Whalen and Reeburgh [1992] 

5 

Panama 

9°N, 80°W 


Keller [1990] 

6 

Canada 

54°N, 105°W 


Valentine (unpublished) 


Table 2: Site-specific model parameters at the 5 test sites 


Test 

Site 

Ro 

v 

v max 

nsoil 

nroot 

T 

veg 

1 

1 

0.60 

45pM*h _1 

80cm 

50cm 

15 

2 

2 

0.30 

20|lM*h' 1 

80cm 


4 

3 

3 

0.34 

20pM*h* 1 

70cm 

30cm 

12 

4 

4 

0.30 

3pM*h _1 

thaw depth 

0cm 

0 

5 

5 

2.80 

20pM*h~ 1 

128cm 

74cm 

9 


Table 3: Sensitivity of simulated CH 4 emissions to climate forcing at site 1 


Run 

1991 

1992 

1993 

T+1°C 

+20.5% 

+20.7% 

+20.5% 

T-1°C 

-16.6% 

-17.0% 

-17.1% 

wt+lOcm 

+25.3% 

+ 1.4% 

±0.0% 

wt-lOcm 

-36.6% 

-20.6% 

-10.8% 

control 

78.2 gCH 4 *yr _1 

70.4 gCH 4 *yr- 1 

65.5 gCH 4 *yr- 1 

















































































